<HTML>
<CENTER><A HREF = "http://lammps.sandia.gov">LAMMPS WWW Site</A> - <A HREF = "Manual.html">LAMMPS Documentation</A> - <A HREF = "Section_commands.html#comm">LAMMPS Commands</A> 
</CENTER>






<HR>

<H3>compute body/local command 
</H3>
<P><B>Syntax:</B>
</P>
<PRE>compute ID group-ID body/local input1 input2 ... 
</PRE>
<UL><LI>ID, group-ID are documented in <A HREF = "compute.html">compute</A> command 

<LI>body/local = style name of this compute command 

<LI>one or more keywords may be appended 

<LI>keyword = <I>type</I> or <I>integer</I> 

<PRE>  <I>type</I> = atom type of the body particle
  <I>integer</I> = 1,2,3,etc = index of fields defined by body style 
</PRE>

</UL>
<P><B>Examples:</B>
</P>
<PRE>compute 1 all body/local type 1 2 3
compute 1 all body/local 3 6 
</PRE>
<P><B>Description:</B>
</P>
<P>Define a computation that calculates properties of individual body
sub-particles.  The number of datums generated, aggregated across all
processors, equals the number of body sub-particles plus the number of
non-body particles in the system, modified by the group parameter as
explained below.  See <A HREF = "Section_howto.html#howto_14">Section_howto 14</A>
of the manual and the <A HREF = "body.html">body</A> doc page for more details on
using body particles.
</P>
<P>The local data stored by this command is generated by looping over all
the atoms.  An atom will only be included if it is in the group.  If
the atom is a body particle, then its N sub-particles will be looped
over, and it will contribute N datums to the count of datums.  If it
is not a body particle, it will contribute 1 datum.
</P>
<P>For both body particles and non-body particles, the <I>type</I> keyword
will store the type of the atom.
</P>
<P>The <I>integer</I> keywords mean different things for body and non-body
particles.  If the atom is not a body particle, only its <I>x</I>, <I>y</I>, <I>z</I>
coordinates can be referenced, using the <I>integer</I> keywords 1,2,3.
Note that this means that if you want to access more fields than this
for body particles, then you cannot include non-body particles in the
group.
</P>
<P>For a body particle, the <I>integer</I> keywords refer to fields calculated
by the body style for each sub-particle.  The body style, as specified
by the <A HREF = "atom_style.html">atom_style body</A>, determines how many fields
exist and what they are.  See the <A HREF = "body.html">body</A> doc page for
details of the different styles.
</P>
<P>Here is an example of how to output body information using the <A HREF = "dump.html">dump
local</A> command with this compute.  If fields 1,2,3 for the
body sub-particles are x,y,z coordinates, then the dump file will be
formatted similar to the output of a <A HREF = "dump.html">dump atom or custom</A>
command.
</P>
<PRE>compute 1 all body/local type 1 2 3
dump 1 all local 1000 tmp.dump index c_1[1] c_1[2] c_1[3] c_1[4] 
</PRE>
<P><B>Output info:</B>
</P>
<P>This compute calculates a local vector or local array depending on the
number of keywords.  The length of the vector or number of rows in the
array is the number of datums as described above.  If a single keyword
is specified, a local vector is produced.  If two or more keywords are
specified, a local array is produced where the number of columns = the
number of keywords.  The vector or array can be accessed by any
command that uses local values from a compute as input.  See <A HREF = "Section_howto.html#howto_15">this
section</A> for an overview of LAMMPS output
options.
</P>
<P>The <A HREF = "units.html">units</A> for output values depend on the body style.
</P>
<P><B>Restrictions:</B> none
</P>
<P><B>Related commands:</B>
</P>
<P><A HREF = "dump.html">dump local</A>
</P>
<P><B>Default:</B> none
</P>
</HTML>
